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Abstract 



The all-orders /3eta function is used to study disordered Dirac fermions in 
2D. The generic strong coupling fixed 'points' of anisotropic current-current 
interactions at large distances are actually isotropic manifolds corresponding 
to subalgebras of the maximal current algebra at short distances. We argue 
that IR fixed point theories are generally current algebra cosets. We illustrate 
this with the simple example of anisotropic su(2), which is the physics of 
Kosterlitz-Thouless transitions. We propose a phase diagram for the Chalker- 
Coddington network model which is in the universality class of the integer 
Quantum Hall transition. One phase is in the universality class of dense 
polymers. 
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I. INTRODUCTION 



Two is the critical dimension for Anderson localization, namely, above 2 dimensions 
states are localized only above a critical strength of the disorder (impurities) whereas below 
2 dimensions states are localized for any strength of the disorder 0,3- For this reason one 
can expect rich phase structures in two dimensions. It is well understood that a complete 
understanding of the Quantum Hall transition involves derealization || [|J . In these local- 
ization problems, understanding the quantum phase transition in the conductivity requires 
finding the right renormalization group fixed point at large distances. This usually occurs 
at strong coupling, which is why these problems are considered difficult. 

The conventional theoretical framework leads to sigma models on various spaces ||. For 
the quantum Hall transition the sigma model was worked out by Pruisken || . These results 
are important, but the sigma models have generally been too difficult to solve. The most 
recent progress can be found in [[7J. 

Critical theories in 2 dimensions are conformally invariant and this imposes strong con- 
straints on the theory ||. Early efforts to use conformal field theory methods in the study 
of disordered Dirac fermions were made by Bernard ||, Mudry et. al. [TIJ, and Nersesyan 



et. al. [ Xlfl - Perturbed conformal field theory methods were also used for disordered sta- 



tistical mechanical models [|T^] generally these are more difficult problems since the 
non-random theories are interacting so that one has to use replicas. Many of the important 
problems here also are driven to strong coupling under renormalization group (RG) flow. 

In this work we describe a general approach to disordered Dirac fermions at strong 
coupling which uses the all-orders /3eta functions proposed in ||14|| . These models are not in- 
tegrable, but one can nevertheless sum up all orders in perturbation theory for the /3eta func- 
tion. In some insightful work, Chalker and Coddington proposed a simple network model 
which is now believed to be in the universality class of the integer quantum Hall transition 



17 1 . This network model has been shown to be equivalent to a certain theory of random 



Dirac fermions and this is the model we analyze using our methods. This model also 



appeared in the work |19| . A sigma model formulation was proposed in |20[ . 

This paper is a first attempt at understanding the implications of these all-orders /3eta 
functions, and for this reason is in part conjectural. The /3eta functions generally do not 
have non-trivial zeros, but rather have poles. We point out that /9eta functions with poles 



are also known to occur in supersymmetric gauge theory [21]. In the next section we propose 



a general scheme for the fixed points of marginal symmetry breaking (anisotropic) current- 
current perturbations for a current algebra Q in 2 dimensions based on these /3eta functions 
and certain hypotheses. Since there is no small parameter in these theories, one generally 
does not flow to a single fixed point, but rather to a fixed point manifold in the couplings. 
In the case of large couplings we show that this manifold generally corresponds to a sub- 
current-algebra 7i of Q and argue that the infra-red fixed point is the current algebra coset 

Q/n. 

The /3eta function shows an interesting duality, namely g flows to 1/g as the length scale 
varies from r to 1/r. This is a new form of duality since the more well-known manifestations 
occur in scale-invariant theories. Here, theories at a scale r and coupling g are equivalent 
to a theory at scale 1/r and coupling 1/g. 

Assuming our hypotheses for interpreting the /Seta functions are correct, they turn out to 
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predict a very tight phase structure; in fact the resulting phase diagram is exact in the sense 
that the phase boundaries are precisely determined. We are unaware of any other models 
where the exact /3eta function allows one to determine so precisely the phase diagram. As an 
illustration we work out the simplest possible case of an anisotropic su(2) perturbation. This 
is the subject of Kosterlitz-Thouless transitions which have previously only been studied at 



weak coupling | 16| . The phase diagram is unexpectedly rich. At small coupling one phase 



exhibits the symmetry restoration recently studied at one loop J^J |£3], however it appears 
this symmetry does not persist at strong coupling. Another phase has a line of fixed points 
corresponding to a free boson at a certain radius of compactification. 

If we extend the physical regime of couplings of the network model, the resulting disor- 
dered Dirac theory again has a number of phases. Though many of the features are similar 
to the anisotropic su(2) case, some features are novel and not completely understood. One 
phase (PSLg) has a line of fixed points related to PSL(1\1) studied in p4| , and appears to 
be the only massless phase. It is in the same universality class as dense polymers and is also 
closely related to the disordered XY model [p5|| . A sigma model version of this theory was 



proposed by Zirnbauer in connection with the quantum Hall transition 41 , based on the 



work j|5| ||46|| . In the physical regime of the network model we show that this phase is not 
easily realized since it corresponds to imaginary gauge potential. The other phases do not 
automatically fall into the category of Q/TC, as explained below. 

An unexpected feature of our analysis is that initially positive couplings which are vari- 
ances of real disordered potentials can flow to negative values, which in turn can be viewed 
as corresponding to imaginary potentials and thus non-hermitian hamiltonians. This could 
signify that our approach breaks down for reasons we do not yet comprehend. If this feature 
turns out to be sensible, it implies that the studies of derealization transitions in non- 
hermitian quantum mechanics made by Hatano and Nelson [27] may have some bearing on 
hermitian systems. 



II. GENERAL STRUCTURE OF FIXED POINTS 

The general class of quantum field theories we consider are current-current perturbations 
of a conformally invariant field theory. We assume the conformal field theory possesses left 
and right conserved currents J a (z) and J a (z) in the usual way [pj ||29|| , where z,z are 
euclidean light-cone coordinates, z = (x + iy)/\^2, z = [x — iy)/\/2. These currents satisfy 
the operator product expansion (OPE): 

J a (z) J\0) = ^-r] ab + -ff J c (0) + .... (2.1) 

z z z 

where k is the level and rj ab the metric on the algebra, and similarly for J. We will refer to 
this current algebra as Qk, where Q is a Lie algebra or superalgebra and the formal action 
for this theory as Sg k . In the superalgebraic case each current J a has a grade [a] = or 1 
corresponding to bosonic verses fermionic. The metric rf b generally cannot be diagonalized 
in the superalgebra case. The metric and structure constants have the following properties: 

rj ab = (— )[ a ][ b ]^ 6a ; f ab = f ba ^ f abc — _(_)[ 6 H C ] f acb {2.2) 

where f abc = f^ b r] lc and i] a bi] bc = 5^. 
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The theory Sg h can be perturbed by marginal operators built out of left-right current 
bilinears: 



s = s 5k + [ E 9 a o\ o A = Yl di,rr (2.3) 

J Z7r A a,a 

The simplest example is a single coupling g with d aa = rj aSi defining the quadratic Casimir 
of Q. For su(N) at level 1, this is the non-abelian Thirring model, or equivalently the chiral 
Gross-Neveu model. We are mainly interested in more general situations where the tensors 
d A b break the symmetry Q, i.e. are anisotropic. 



In |T4| an all-orders /3eta function in a certain prescription was proposed and we now 
summarize this result. The theory is not renormalizable for any choice of d A . The three 
conditions ensuring renormalizability are: 

(-) m K b d B c jrf? = C£ B d% (2.4) 

= D c B d° ba (2.5) 

dfjftft = R A b V ac df b (2.6) 
The first condition is equivalent to closure of the operator algebra of O a : 

O a (z, z)O B (0) O c (0) (2.7) 

zz c 

and guarantees one-loop renormalizability. The other two conditions are necessary at two 
loops and higher. The structure constants D, R are also related to an OPE. Define the 
operator T A built out of left-moving currents only and d A : 

T A (z) = d A J a (z)J b (z) (2.8) 

The operator T A is a kind of stress tensor; in the isotropic case it is the affme-Sugawara 
stress tensor up to a normalization. One finds 

T A (z)O B (0) ~ (2kD AB + R^D BD ) O c (0) (2.9) 

Given a particular theory, the above equation is an efficient way to compute R, D; the D 
term is distinguished from the RD term by being proportional to k. The conditions required 
for T A to satisfy the OPE of a consistent stress-tensor, the so-called master equation studied 



extensively in p!5| , involves the same objects C, D, R, but appears to be a stronger than our 
renormalizability conditions. 

The renormalization group (RG) structure constants have the following properties: 

D AB = D BA , C AB = C BA } D AC D D E B = D AB D D E C (2.10) 

which can be proven from the defining relations or from the OPE's. 

Let us arrange the couplings into a row vector g = (gi, g 2 , ...). Since the higher loop 
expansion is in kg, let us define 

9 = kg/2 (2.11) 
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Let D(g) be the matrix of couplings 



D{9)b = Y. D b9c (2.12) 

c 



Given two row vectors v ' we define a new row vector from C: 



c{v'y) A = Y,vW c c B A c (2.i3) 

B,C 



Finally, define 



' g '-' g T^DW (2 - M) 

Then the /3eta function can be expressed as 

/% = = I {-\c{g' ,g'){\ + D 2 ) + C(g'D,g'D)D - g'DRD^j (2.15) 

where r is a length scale and D = D(g). The flow to the IR corresponds to increasing r. 
For the purpose of studying fixed points, it will be convenient to define 

C£ B = R^)Dq D (2.16) 

Note that C appears directly in the OPE ( |2.9|) and the /3eta function so there is no need to 
compute R by itself. One can show 

C(v\v 2 D(g))=C(v\v 2 )D(g) (2.17) 

using fl2.10|) . The /5eta function can then be arranged into the form: 

Pa=\[- C @M l ~ D ) 2 + & W+,g-) - C(g„,g + ) (2.18) 
- ((C + C)(s-,g + )) D-((C + C)(g + ,g_)) d] (2.19) 

where again D = D(g) and we have defined 

& =»T±-b) (2 ' 20) 

We now study the possible fixed points based on the above /3eta function. In the specific 
examples we have studied has no non-trivial zeros for any finite values of g. The only 
zeros of /3j in the two examples below are trivial in the sense that they correspond to g = or 
exactly marginal directions, such as u(l) directions. We believe this is the generic situation 
but cannot rule out other examples with non-trivial zeros. 

The /?eta function generally has poles as can be seen from ( 2.14 ). These arise from 
summing geometrical series in perturbation theory. Since these series do not converge for 
g > 1 one should question their validity for large g. This being the first attempt to make 
sense of these /3eta functions, in the sequel we will simply assume the above /5eta function 
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is valid for all g, which amounts to an analytic continuation of the perturbation series. As 
we will see, the resulting RG flows can be interpreted at large g. 

Assuming the /3eta has no non-trivial zeros, then this implies that we must look for fixed 
point manifolds. Suppose all couplings g have large absolute value. Since g' ~ 1/g, g± ~ 1 
and D(g) ~ g, the /5eta function is dominated by the C + C terms which are of order g. For 
large g we therefore expect the possibility of flowing to manifolds in the coupling constant 
space satisfying 

(C + C)(g„,g + ) = (C + C)(g + ,g„) = (2.21) 

In the examples below the solutions to the above equations are the same as for 

C(g,g)+C(g,g) = (2.22) 

The above argument is not strong enough to conclude that the only asymptotic flows are 
solutions to ( |2.22| ). Rather, we expect that solutions to ( |2.22j ) represent a generic class of 
fixed points. This is strongly supported by the numerical analysis of the exact /3eta functions 
in the examples below; it is found that when couplings flow to infinity, one indeed generally 



flows to solutions of the equation Q2.22 



Solutions to the above fixed point manifold equation are generally couplings that corre- 
spond to invariant subalgebras of Q. Let Q a denote generators of Q: 

[Q a ,Q b } = ffQ c (2.23) 
In the conformal field theory these are realized as conserved charges: 

Suppose the d A b define invariants of Q: 

Q a ,d A Q b Q c ]=0 (2.25) 

Then this implies 

dfjf = (-) [a][c] dif b ai (2.26) 



Using this in the renormalizability conditions ( |2.4| ) one finds that this implies C + C = 0. 

Interesting possibilities for non-trivial fixed points thus arise when the d A couple com- 
muting subalgebras of Q. Let 

A = VabJPl (2-27) 

where J a are the currents for a subalgebra Ti A of Q, r\ A defines the Casimir for 7i A , and these 
subalgebras commute for A ^ B. The currents J a satisfy a current algebra TC A , where the 
level fvA is generally not equal to k. This is sometimes referred to as a conformal embedding. 
Only the currents J a matter in ( |2.25| ) and the fixed point manifold condition C + C = is 
satisfied. 

Let us further suppose that each 7i A has a unique quadratic Casimir so that 



6 



Vtjftfil = C^5 c d (2.28) 

where is the quadratic Casimir in the adjoint representation of TL A . Not all algebras 
have this property; in particular the superalgebras with indecomposable representations 
which we will encounter below do not. Assuming Q2.2S ), one finds that 



C A A - ~ C adji R B - C adj$B, D A A = 1 (2.29) 

The /3eta functions then obviously decouple and take the simple form: 

- ^orb (230) 

The /Jeta function ( |2.30| ) has some interesting properties. Let us focus on a single coupling 
qa = g- Define a dual coupling 



9* = i/g (2-31) 



The /9eta function satisfies 



P*(9*) = -(3(9^9*) (2.32) 

This implies that 

Sir) = (2.33) 

This behavior can be seen from the solution: 

(g - l/g)/2 + hgg = ^- log r/r (2.34) 

where ro is the scale where g — 1. As shown below, the anisotropic case also exhibits this 
kind of duality. 

There are four cases to consider. First suppose C a dj/k > 0. Then g always increases. 
From the duality ( |2.33| ) one sees that g > at short distances (small r) flows to g = oo at 



large r. On the other hand any g < flows to 0~. Thus g > is marginally relevant and 
g < marginally irrelevant. This confirms the expectations based on one-loop. Note that 
the double pole in the /3eta function at g = — 1 has no physical effect whatsoever: near the 
pole g grows faster as r Cad ^ 2k rather than logr but simply flows through the pole. When 
Cadj/k < everything is reversed: g > flows to + and g < flows to — oo. Thus in this 
case of isotropic couplings, the /9eta function appears to be physically sensible for all g, and 
this supports the validity of the analytic continuation of the /3eta function. 

We thus propose the following generic scheme for non-trivial infrared (IR) fixed points 
that can occur when couplings flow to infinity. For every marginally relevant g^, gA flows 
to either ±oo. The degrees of freedom coupled by gA are thus massive and should decouple 
in the flow to the IR. This picture was used in the work [OT on the spin quantum Hall 
effect. For example consider an su(3)i current algebra perturbed by the marginal operator 
that isotropically only couples an su(2)i sub-current algebra with coupling g. The sw(3)i 
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current algebra can be bosonized in terms of two free bosons 0i,02 and has c = 2. The 
su(2)i current algebra can be bosonized in terms of a single scalar field which is a linear 
combination of 0i, 02, and has c = 1. As g goes to oo under RG flow, correlation functions of 
tend to zero. Thus the field decouples and the IR fixed point is su(3)i/ su{2)\. Generally, 
the fixed point is thus 



n 

IR fixed point = U k = ® A Ut A (2.35) 



where <S>a is only over marginally relevant The IR theory is non-trivial only if H, ^ 
Q. This situation corresponds to a massless phase and Qk/'Hk are the massless degrees 
of freedom. When Qk/'Hk is empty, e.g. when Q = TC, then this is a massive phase. 
In localization problems a massive phase corresponds to all states being localized. Coset 
conformal field theories were studied in generality in |31|. The Virasoro central charge 



at the fixed point is (9 — c n , and the conformal scaling dimension A of an operator at the 
fixed point is 

A g/n = A Q _ A n ( 2 36 ) 

An interesting open question is whether other models can flow to the more general solutions 
of the master equation |0| which do not correspond to cosets. 



Due to the existence of poles in the /3eta function, not all flows that we find in the 
examples below fall neatly into the above scheme. In some cases below the RG flows are 
attracted to these poles. In most of these cases, as one approaches the pole, some of the 
other couplings flow off to infinity in such a way that suggests one can still perform a coset. 



III. ANISOTROPIC SU(2) 



In this section we illustrate the scheme of the last section in the simplest possible ex- 
ample of an anisotropic su(2). This is the subject of Kosterlitz-Thouless flows which have 
previously only been studied at weak coupling [nj. We originally worked out this example 
to check of the validity of the exact /3eta functions. The phase diagram is unexpectedly rich. 

We normalize the currents as follows: 



and consider the action 

S = Ssu(2) 

A simple computation using 



±-^(0), 

z 



J + (z)J-(0) 



2791) gives the RG data: 



C{ 2 
Dl 1 



Cf 

ly 2 



-2 



Cf 



1/2 



2z l z 



(3.1) 



(3.2) 



(3.3) 
(3.4) 
(3.5) 



The matrix D(g) is diagonal: 



S 



D{9)=[H 2 °,o) ( 3 - 6 ) 



g 2 /2 



The resulting /3eta functions are 



P91 (1-A;V/16)(1 + W4) l ' J 



Pa2 (1-A;V/16) 2 



There is again an interesting duality in these /3eta functions. Define the dual couplings 
Then the /3eta function satisfies 

W) = ^ <?*) (3.10) 

The above /3eta functions predict a tight but unexpectedly rich phase diagram. The 
duality (|3.10| ) explains some features of the diagram, namely, the self-dual lines g = g* are 
phase boundaries or lines of attraction. The solutions to C + C = are (i) gi = g 2 , which 
corresponds to the subalgebra Tik = su{2)k and (ii) g± — corresponding to TC — u(l). 

The structure of the phase diagram is determined in part by the behavior near the poles 
at gi,g 2 = ±4/&. Consider the pole at (gi,g 2 ) = (4, 4)//c and let gi j2 = 4(1 + e\ j2 )/k. Near 
e = the behavior is 

A.-S-^, /5 e2 ^^| (3.11) 
e i e i 

Around this pole, this leads to the behavior shown in figure 1. In this figure, heavy lines 
correspond to phase boundaries. The region, or phase, A is attracted to the line g\ = g 2 , 
since this line is stable in region A. Beyond (g±,g 2 ) = (4,4)//c the line g\ = g 2 becomes 
unstable. Since the line becomes unstable beyond the pole, one has to reach the line exactly 
in the region A before flowing off to infinity along it; otherwise one can flow elsewhere (see 
below). If one is on the line, then one flows to infinity and the IR fixed point is the empty 
coset su{2) k / su{2) k . 

In region D, g 2 flows to oo whereas g\ flows to a constant, but after a finite scale 
transformation. Because one coupling is blowing up, the flow cannot be continued to larger 
scales numerically. Since the ratio g\j g 2 flows to zero, one possible interpretation of this flow 
is that gi is effectively zero and g 2 = oo. The subalgebra coupled by g 2 is H = u(l), and the 
IR fixed point would thus be su(2)k/u(l). This is a well-known conformal embedding [|32] and 
corresponds to the critical theory of Z k parafermions with central charge c = 2(k — l)/(k+2), 
where k = 2 is the Ising model. However we emphasize that since the RG flow cannot be 
continued to arbitrarily large length scales, this is not a true fixed point. 

In phase E, g\ grows to infinity and g 2 to a non-universal value < g 2 < 4. Here the 
flow can be continued to arbitarily large scales. For fc = lwe interpret this as a sine-Gordon 
phase. At k = 1 we can bosonize the currents 
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where ip(z) is the ^-dependent part of a free massless scalar field 0. Viewing the g 2 coupling 
as a perturbation of the kinetic term and rescaling the field cb one obtains the sine-Gordon 
action 



S = — I d 2 x 

47T 



(3.13) 



The region A is also sine-Gordon like at small coupling. The operator cos (3 has scaling 
dimension T = (3 2 . To relate (3 to g 2 , consider the /3eta function when g± is small. Then 
Pgi = 9i92/ (1 + ^2/4) • Since this is proportional to gi, we can identify the dimension of the 
coupling as dim(<7x) — 52/(1 + = 2 — T. Thus: 

In general, we expect (3 to be a function of both gi, g 2 . Note that for small g±, values of g 
between and 4 correspond to < p 1 < 2 which is the expected regime for the sine-Gordon 
model. (The conventional sine-Gordon coupling (3 is (3 = f3y4ar.) For k > 1 the appropriate 
generalization is the fractional super sine-Gordon model |34[ . 



In the G phase g\ flows to zero and g<i to a fixed value on the line —4 < kg 2 < 0. For 
example, as r — > 00 one finds 

((?!, 2 ) = (2, -3) — > (0, -2.337693444..) (3.15) 

when k = 1. This line of fixed points corresponds to a free boson at some radius of com- 
pactification determined by the value of g 2 on the fixed line. 

In region I one is attracted to the isotropic line and then flows toward the poles at 
(%x,%2) = (4,-4). 
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FIGURES 




FIG. 1. Phase Diagram for anisotropic su(2). Phase boundaries are solid heavy lines. Heavy 
dashed lines denote lines of fixed points. Shaded circles are fixed points. 

The flows toward the poles in region A and I is rather delicate, and the outcome cannot 
be studied numerically by solving the differential equations because of the singularity. As 
stated above, one possibility is to continue to flow off to infinity along the isotropic lines 
once one reaches the poles. However this appears to be inconsistent with the duality ( p.lO ) . 
The duality implies that if a flow passes through the self-dual point g = g*, then g and g* 
are on the same RG trajectory, where g flows to g* as r goes to 1/r. We can use this to 
extend the flows through the poles^. Namely, the region in A below the isotropic line flows 
to the region D, whereas the region above the line flows to region E. Similarly, the region 
above the isotropic line in I flows to region G, whereas the region below the line flows to A 
then to D. This implies that the symmetry restoration seen in region A is destroyed once 
the flow reaches beyond the pole. 

The remaining phases are a mirror image of the above. Though the line g\ = — #2 does 
not appear as a solution to C + C = 0, and it would naively break the su(2) symmetry, 
it turns out that it does possess an su(2) symmetry. There is an automorphism of su(2) 
J ± — > —J^ that preserves the algebra. We can perform this automorphism on the left 
currents only. This takes O 1 — > — O 1 , i.e. flips the sign of g\. The phase diagram then 
has symmetry with respect to reflections about the gi axis. Thus the C, F, B, I phases are 
similar to the D, E, A, H phases. 



An earlier version of this paper did not consider this possibility, which will be further reported 
on in |p3[ . 
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A check of the above /3eta function is based on comparison with the sine-Gordon 
/3eta function which is known to two loops. For a recent discussion, see [35]. These per- 
turbative computations were performed around (3 = a/2 and g% equal to zero. Define new 
couplings 6, a as 



(3 2 1 - g 2 /4 
2 1 + #2/4 



(3.16) 

Then the /3eta function (|3.7|) implies the following /3eta functions for 5, a to 5 loops: 

n ^ 2 1 Or 1 4 1 2 r 9 ^ 4r 

da = a a o a a o a o 

H 32 16 4096 32 2048 

p a = -2a5 - ^-a 3 - ^-a 3 5 - — \— a 5 (3.17) 



It is known that with more than one coupling, not all 2 loop contributions are prescription 
independent. Define new couplings to second order as 



5' = 5 - —a 2 
32 



a = a 



(3.18) 



Then the /3eta function (|3.17|) to two loops leads to 

Ps> = 



1 /2 i 1 12 cl 

— a H a o 

32 16 



(3 a , = -2a'5' - -a' 3 
64 



(3.19) 



This agrees with the result presented in |3(J [35|. 



IV. THE NETWORK MODEL 
A. /3eta function 



The network model can be mapped onto a model of disordered Dirac fermions JTS) . The 
two-component hamiltonian in the two spacial dimensions x, y is 

H = {-id x - A x ) <r* + ^7= (~idy -A y )a y + V + Ma z (4.1) 

where a are Pauli matrices and A(x, y), V(x, y), M(x, y) are real disordered potentials. Dis- 
order in A, V and M respectively corresponds to randomness in the individual link phases, 
the total Aharanov-Bohm phase per plaquette, and the tunneling at the nodes. This model 
was also considered in the work ||19|| . 

Introducing 2-component Dirac fermions 



*=[f), ** = (V_,V-) (4-2) 



one needs to study the action 
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s 



d 2 x 



d 2 x 



tfj_(d z - iA z )ip + + ip-(d B - iA z )ip + 



(4.3) 

(4.4) 
(4.5) 



where z = (x + ij/) /y/2, z = (x - iy) fy/2 and A z = (A x - iA y ) /y/2, A z = (A x + %A y ) /y/2. 
We take the potentials to have the following gaussian distributions: 



1 



d 2 x. 



P[M] = exp 
P[A] = exp 



2^, 



d 2 x 



2ix 



-M 



1 rd 2 x 



9a 



2tt 



A z A z 



(4.6) 
(4.7) 
(4.8) 



The couplings g v , g m , g a are positive variances of the potentials. If A — * iA, then g a — ► — g a , 
and similarly for M, V. Thus negative g's can be interpreted as corresponding to imaginary 
potentials. 

Since we are dealing with a free theory, we can use the supersymmetric method for dis- 
order averaging |3"7fl . In the present context this method was studied in || jlTJ. Introducing 
bosonic ghost partners (3±, /3± of the fermions and performing the gaussian integrals one 
obtains: 



S eS = S bee + / ^ (9vO v + g m O m + g a O a ) 



S'free is a free c = conformal field theory with the action 



S 



free 



d 2 x 



d z ip+ + + P-9 Z P+ + P-d g p+) 



(4.9) 



(4.10) 



A treatment of the c = — 1 ghost system can be found in [|38f 



The maximal conserved currents of Sf ree are the 8 possible bilinears in the fermions and 
ghosts: 



J = P+P- > J± = PL S ± = ±^ ± (3-_ 



Using the OPE's: 



^ + (z)V>_(0) ~ ^-(^#+(0) ~ 1/z 
P + (z)P-(0)~-/3-{z)P + {0)~l/z 

one finds that the currents satisfy the osp(2\2)k=i current algebra: 



S± = ip T P T 



(4.11) 



(4.12) 



J(*)J(0) ~ ~ 



H(z)H(0) ~ A 
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J(z)J±(0) ~ ±- J±, 
J(z)S±(0) ~ ±^S ± , 
H(z)S+(0) ~ ±-£±, 
J±(z)S T (0) ~ ~£±, 



J + ( 2 )J_(0)~^--J 
J(z)£±(0) ~ ±-£± 

#(42±(o) ~ t-£± 



J±(«)^(0) 



(4.13) 



1 



S ± ( 2 )5 ± (0)~±-4 

z 

S + (z)S-(0)~\ + l(H-J) 

z z z 

§ + (z)S-.(0)~~ + ±(H + J) 
z z z 

To better reveal the algebraic structure, we define new couplings 

g v O v + g m O m = g + + + g-O- 

with 

g ± = g v ±g m , ± = (O v ±O m )/2 

The perturbing operators can then be written in terms of currents in the following 

1 



- ( / / ' ' 



+ = S + S- - 

0~ = JJ -HH + 3+3- - SS 
O a = (J-H)(J-H) 



Using the OPE's ( |4.13|) , it is straightforward to compute the RG data from 
One finds the non-zero values: 



and 



a+ H 



C++ = c, 



(4.14) 

(4.15) 
way: 

(4.16) 

3 ID- 

(4.17) 



C++ 



-c: 



c: 



-ci + 



-c+ 

= -4 



-2 



(4.18) 



The non-zero -D's are 



D a - = D- a = DZ~ = -D+ + = -1 
The matrix D(g) is thus non-diagonal. In a basis g = (g+, g_, g a ): 

D(g) 



(9+ 











-9- 




vO 





-9- 
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(4.19) 



(4.20) 



After some algebra we find the /5eta functions: 

8g + {g\{2g a -g-+2) + 2g_(2 - g_) + 8g a 



(A-g 2 + )(2-g_y 

{±-9. 

4 Ug\ - «£)(16 - glgl ) + ±g a g\ (2 + y_)(2 - S-) 2 



= 7/ I 2 V (4-21) 



(4-^)2(2-^)2 



B. Phases of the network model 



The /3eta functions (|4.21| ) have a precise phase structure that is not readily apparent from 



their form. We studied the phase structure by analyzing the behavior in the vicinity of the 
poles combined with some modest numerical work. Namely, the /3eta function differential 
equations were solved numerically for g(r) for a variety of points in each phase. 

To begin describing the phase diagram we start with the solutions to C + C = 0. There 
are two solutions: (i) g + = —g~,g a = and (ii) g + = 0. The first solution corresponds to 
g v = 0,g± = ±g m and the perturbation is thus g m O m with 

O m = (-JJ + HH- S + S_ + S_S + + S+S- - S-S + + \{J+J- + J-J+)) (4-22) 

O m is built on the Casimir of osp(2\2), so this corresponds to an osp(2\2) symmetric manifold. 
The /3et& function is 

^ <0h (423) 

When g a ^ this line is unstable so it does not play a significant role in the phase diagram. 
However as we will see, there is one regime that may be attracted to it. 

For the other solution g + = with g~,g a arbitrary, the only currents in the perturbation 
are J,H,S± which generate gl(l\l) at level k — 1. The two operators 0~ , O a correspond 
to the two independent Casimirs of grZ(l|l), reflecting the indecomposability of the adjoint 
representation. The /3eta functions reduce to 

= - (2- g l_r Pb ~ = ° (4 ' 24) 



This model was studied in |^| in connection with the disordered XY model and the localiza- 



tion problem of electrons randomly hopping on a lattice with ir flux per plaquette . The 



coupling g + = corresponds to g v + g m — which means that one of the couplings g V)m is 
negative. In [24"] such negative couplings arose naturally in a different hermitian hamiltonian 



with twice as many degrees of freedom. In the model we are studying, as we will see, a large 
regime of couplings, including initially all positive couplings, are attracted to this 5^(1 11) 
invariant manifold. This is rather unexpected: starting from a hermitian hamiltonian with 
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positive variances g, in some regimes the couplings flow to negative values which correspond 
to imaginary potentials and non-hermitian hamiltonians. 

There exists another osp(2\2) invariant line which does not appear as a solution to 
C + C = 0. It arises due to the automorphism of osp{2\2): 

J± -> -J±, S± -> -S± (4.25) 

which does not change the algebra. Performing this automorphism on the left-moving cur- 
rents only sends O m — ► — O v . Therefore the line g + = g~,g a = 0, corresponding to g m = 0, 
is also osp(2\2) preserving with the /3eta function: 

= i& (4 - 26) 

For g v > 0, this is marginally relevant and is thus a massive theory. It's an integrable theory 
and the exact S-matrix was proposed in [39|. In the network model, one phase (called O 



below) appears to be attracted to this line but with g v < 0, and thus flows to zero, and is 
massless. 

As in the su{2) case the global features of the phase diagram largely, but not completely, 
follow from the behavior near the poles g+,g~ = ±2. Consider first the vicinity of the pole 
(g + , g_) = (2, 2). Letting g± = 2 + e±, one has 

8 8 «2H 3 (4 27) 

Pe+ e + e 2 _ ' ^" 4' P9a 4 16e 2 1 J 

If ^ a > 0, then e + is attracted to zero, whereas e_ always grows. Here, g + flows to 2 and 
^_ flow to infinity if g_ > 2. Since g+/g~ flows to zero, we interpret this phase as flowing 

along the g/(l|l) invariant line. 

The behavior around (g + ,g_) = (±2, —2) is different. Here one finds: 

&_«^, (4-28) 
e + e + e + 

Again when ^ a > one is attracted to e + = 0. However when e_ = 0, 8 e _ = 0. 

A simplifying feature is that there are no poles in g a . For each phase, given initial values 
for g±, we find that for some value g a0 which depends non-universally on g±, then g a > g a0 
and g a < g a0 are distinct phases. Examining (3 g+ near the poles g + = ±2, one finds 

9ao « (gl - 4)/8 near g + = ±2 (4.29) 

For large g + , (3 g+ w -8g + (2g a - g-)/g 2 __, thus 

g a0 ng_/2 for # + >l (4.30) 



We list the regions with distinct behavior below. The phase diagrams are shown in 
figures 2,3. For each phase we determine the density of states exponent. To study the 
average density of states p(E) we shift H — > H — E leading to a coupling E$e in the 
effective action with 
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(4.31) 



The density of states is proportional to the one-point function of §e'- 

pjE)<x(<l> E ) (4.32) 

Let T E denote the scaling dimension of the operator in the IR. Then since E has scaling 
dimension 2 — T E , 

~p(E) oc E Te ^ 2 - Te \ asE^O (4.33) 



g a > gaophases 

(i) dXY ± phases. Here g + flows to ±2 after a finite scale transformation, and g>_ , g a flow to 
+oo. This is similar to the D region of anisotropic su(2). The ratio g+/g~ — > 0, so consid- 
ering g+ as effectively zero, one interpretation of the flow is to the coset osp(2|2)i/p/(l|l)i. 
As in the su(2) case this is not a true fixed point since the flow cannot be extended to 
arbitrarily large length scales. The stress tensor for the g/(l|l)i current algebra conformal 
field theory has a structure that parallels the structure of 0~ , O a . Namely, the stress tensor 
is an affine-Sugawara construction built on the sum of the two independent Casimirs [[5J 

T g im k = -^ j2 - R2 + s + s - - s - s +) + ^ J - H ? ( 4 - 34 ) 

This stress tensor has c = so the coset osp{2\2)\/ gl{l\l)\ also has c = 0. The fermions ip± 
have conformal dimension A = 1/2 with respect to the gZ(l|l) thus the gl(l\l) dimension 
of is 1. This implies = and a constant density of states near E = 0. Though 
the algebra gl(l\l) is smaller than osp(2\2), the ^Z(l|l) dimensions of fields are the same as 
for osp(2\2)i. One can in fact construct a level-1 representation for both of these current 
algebras using the same number of fields: two bosons and a complex fermionic scalar. (See 
e.g. [0].) Thus osp{2\2)i/ gl{l\\)i is empty and this would be a massive phase. 

(ii) ± phases. This phase is characterized by {g +1 gJ) flowing to the pole (±2, —2) while 
g a flows to zero. Unlike previous examples, here one is not attracted to the isotropic lines 
g + = ±g_ before reaching the pole since the line is unstable. The fate of the flow once it 
reaches the pole is again delicate. Based on the example of su(2) it seems most likely that 
this flow spills into the Q ± phases. However here we do not have the duality arguments to 
support this. 

(iii) Q phases. In this region one flows to the poles at g + = ±2 after a finite RG time 
but unlike the dXY ± phases, numerical integration indicates that the other couplings g~,g a 
do not flow to infinity but rather to some finite non-universal values as one approaches the 
pole. Since none of the couplings are flowing to infinity, we cannot interpret this as a coset. 
This phase clearly requires further investigation. 
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FIG. 2. Phase Diagram of the network model with g a > g a Q. Phase boundaries are solid heavy 
lines. Shaded circles are fixed points. The cone opening to the right is the physical regime of the 
network model. 



ga < gaophases 

(i) PSL g phase. Here g + — > 0, but g_ flows to a fixed, non-universal value. g a flows very 
slowly to — oo. In this phase g a o ~ 0. Unlike the flows to the poles, in this case the flows 
can be run to arbitarily large length scales. For example 

(g+, g-,g a ) = (1, 1, -1)— (0, 1.1508434.., -oo) (4.35) 

as r — > oo. In this phase we are again on the gl(l\l) invariant manifold but in contrast 
to the dXY phase there is a line of fixed points corresponding to the value of g>_. This is 
consistent with the /3eta function (|4.24| ) since (3 g _ = 0. 
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gx- 



gx + 



g+ 



§a <g a0 

FIG. 3. Phase Diagram of the network model with g a < g a $. Phase boundaries are solid heavy 
lines. The heavy dashed line is a line of fixed points. 

The massive decoupled subalgebra TC is then generated by only H—J. However it appears 
inconsistent to divide only by H — J. The current H — J is not primary with respect to the 
stress tensor T = (J — H) 2 /A, i.e. T(z)(H — J)(0) ~ 0. A consistent conformal embedding 
is based on the subalgebra Ti = tt(l) <8> u(l) with stress-tensor 



- (h 2 - J' 



(4.36) 



The fixed point is then: 



PSL„ phase : IR fixed point 



osp(2|2)i 

u(i)®u(i) 



(4.37) 



with c = —2. This is a massless phase, since the above coset is not empty. In the formal 
limit g_ — » oo the above fixed point is the same as for the dXY phases. As described in 
2M, dividing by u(l) <S> u(l) leaves the free field theory of a complex fermionic scalar, with 



action 



S 



X 



(4.38) 



This theory is equivalent to a (£,//) fermionic system of conformal scaling dimension (0, 1) 
||38||. It is perhaps the simplest logarithmic conformal field theory |42j| . Interestingly the 
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latter theory was used to describe the dense phase of polymers by Saleur [ji3|j . This theory 
is also closely related to the PSL(1\1) sigma model introduced into the context of quantum 
Hall transitions by Zirnbauer fHfl . The coupling g_ naively corresponds to the exactly 
marginal direction 5S = g^S where S is given in ( f4.38|) , so that g_ is similar to a radius 
of compactification for a free boson. However the modular invariant partition functions 
at c = —2 do not have a continuous parameter, but rather are related to Coulomb gas 
partition functions at certain discrete radii |44||. This suggests that as g a — > oo, g_ can only 



take discrete values in order to lead to a modular invariant partition function. The spectrum 
of anomalous dimensions is then the same as in In particular the twist fields of this 
theory have fixed dimension —1/4 and 3/4 [|2]. 

(ii) gX ± phase. Here g + — ► oo, g a —* — oo, and g_ flows to a finite non-universal constant 
as r goes to oo. For example: 

{9+, 9-, 9a) = (10, 5, 1)— ►(«>, 5.8829077.., -oo) (4.39) 

Though g a flows to — oo, the ratio g a /g+ flows to zero. Since g~/g+ also flows to zero, we 
interpret this as g a = g- =0. 

The coupling that goes to infinity, g + , couples the currents S±,J±. These do not form a 
closed subalgebra, and for this reason this phase is more difficult to comprehend and requires 
some speculation. It could simply be a massive phase corresponding to osp(2\2)i/osp(2\2)i 
since the above currents close on the whole of osp(2\2), as in the massive sine-Gordon phase 
of su{2) (region E). Alternatively let us suppose that only some of the currents are set to 
zero by g + going to oo. Since the commutator of J± with S± closes on gl(l\l) currents S±, 
and the gl(l\l) coupling g_ is not flowing to infinity, we cannot consistently set both S± and 
J± to zero. To distinguish S and J± one may need some further osp(2\2) symmetry breaking; 
this could come from the fact that the operator <3>e breaks osp(2\2). Two subalgebras of 
osp{2\2) involving the above currents correspond to another g/(l|l) generated by (S±, J, H) 
and su(2) ® generated by (J±,J,H). Let us suppose the g/(l|l) is set to zero by 
g + — > oo. As for the dXY phase, the coset osp(2\2)i/ gl(l\l) 1 is empty and this would be a 
massive phase. 

Let us consider the other possibility. The currents J, J± generate an su{2) at level 
k = —1/2. To see this, let J -> 2 J, J± -> ±2 V / 2J ± . Then the new currents satisfy the OPE 
( jOf ) with k = —1/2. The coset in this case is osp{2\2) / su{2) u{l) . The central charge 
of su{2)k is c = 3k/(k + 2) which gives c = — 1 for k = —1/2. Since c = 1 for the u(l), the 
above coset has c = 0. The left-moving conformal dimension of primary fields of spin j in 
su{2) k current algebra is 

A (i) = JU + 1) 
k k + 2 

which equals 1/2,4/3,5/2,... for k = —1/2. The ghost fields (3± transform in the spin 1/2 
representation. Thus in the IR the operator (3 + (3_ has dimension A osp ( 2 ' 2 ) — A su ( 2 ) = 0. 
Dividing by u(l) also leads to dimension zero for the fermion part of Thus again 

Te = and p(E) is a constant. Since the su(2)_i/ 2 is built directly from the ghost fields 
with the ghost fields in the spin 1/2 representation, the osp(2|2)i theory appears to be 
equivalent to sit(2)_i/2 <B> u(l) if the only primary field of su(2)-yz is j = 1/2. The coset 
osp(2\2)i/ su(2)-x/2 <S> u(l) is thus empty and this is also a massive phase. 
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C. Physical regime of the network model 



The physical regime of the network model is g v ,g m ,g a all positive. In the (g + ,g-) plane 
this regime of couplings corresponds to the intersection of g + > g>_ and g + > —g~, which is 
the 90° cone symmetric about the g + axis. For reasons we do not yet understand, couplings 
can flow in or out of this cone. Let us suppose that the initial couplings are in the cone, and 
furthermore that initially g a > 0. The phases dXY + , Q + and + are easily accessible since 
g a is positive and can easily be chosen greater than g a Q. 

Consider now the g a < g a o phases. For PSL g , g a0 is (g 2 _ — 4)/8 near g + = 2. But since 
g 2 _ < 4 in the cone, this phase appears inaccessible to the network model. If a negative g a 
turns out to be physically sensible, perhaps for the reasons described in [fT7[] , then this phase 
can be realized. On the other hand, g a < (g^_ — 4) /8 is easily satisfied with a positive g a in 
the cone for the gX + phase. An example is equation ( [4.391 ). 

In summary, the physical regime of the network model can flow to one of 4 different 
phases dXY + ,0 + , gX + , and Q + . Initially weak couplings flow to the phase Q + , which, 
as stated above, cannot be interpreted with our hypotheses since none of the couplings are 
flowing to infinity. 

In the above phases dXY + and gX + , the density of states is not critical, i.e. p(E) oc E° 
since Te = 0- This is in accordance with the expectation that the disordered 1-copy theory 
is not critical. The conventional wisdom is that one needs to study the 2-copy theory and 
compute disorder averages of the product of retarded and advanced Green functions. This 
can be done using the methods of this paper. Since the 1-copy theory is contained in the 
N-copy version, one expects on physical grounds that the /3eta functions are the same for 
all N. We checked that this turns out to be the case and is a consequence of the zero 
super- dimension of osp(2N\2N). Though the /?eta functions are the same as we described 
for N — 1, the fixed points are different since the current algebras involved are different, 
in particular, osp(2\2)i is replaced by osp(2N\2N)i. This will be described in a separate 
publication. 

Though the localization length exponent v is currently beyond our understanding, let us 
examine the possible exponents for the dense polymer phase (PSL g ). Tuning through the 
critical point generally corresponds to a perturbation 



with e ~ for some operator $ e . Let r e denote the scaling dimension of $ e at the infra-red 
fixed point. The mass dimension of e is then 2 — r e . As e — > 0, there is thus a diverging 
length scale 



The fields in the dense polymer theory all have dimension n/16 plus an integer. In particular 
there is a field of dimension T = 25/16, corresponding the conformal dimension h = 25/32, 
which is a descendent (1+) of the field which is the dense polymer 1-leg operator (h = —3/32) 
times a twist field (h = —1/8) from the x-theory sector of another copy. Note that this field 
does not exist in the 1-copy theory. Such a field leads to v — 16/7. Within this scheme, this 
appears to be the closest one can come to the numerical value 2.35 ± .03 pE9fl . Unfortunately 
we have no further arguments supporting the significance of this operator. 




i oc e 



!/ = l/(2-r c ) 



(4.40) 
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V. CONCLUSION 



Under the hypotheses outlined above we have interpreted most of the RG flows based 
on the all-orders /5eta functions proposed in ||14|| . Much remains to be further clarified, 
in particular the flows that are attracted to poles in the /3eta function after a finite scale 
transformation need further investigation. In any case, this work shows that disordered 
fermions in 2D at strong coupling can have a rich phase structure. In the physical regime 
of the network model the important phases that we could understand as cosets appear to 
have a constant density of states, in accordance with the conventional understanding. In 
addition, our analysis of anisotropic su(2) could have implications for Kosterlitz-Thouless 
transition physics at strong coupling. 

The only certainly massless phase we found for the network model is the c = —2 conformal 
field theory osp(2\2)i/u(l) (g) u(l), which is known to correspond to dense polymers. There 
exists a classical percolation picture for the quantum Hall transition [ i8| , and it is known 



that percolation and dilute polymers are closely related c = conformal field theories |13| 
Our work seems to suggest that with strong disorder the classical dilute polymer theory 
flows to a dense polymer phase. 

We considered the simplest case of one copy of Dirac fermion. In the theory of disorder, 
for the same model one needs to study more copies to compute averages of products of 
correlation functions, and this generally has multifractal behavior. For N copies the short 
distance unperturbed theory has osp(2N\2N)i current algebra symmetry. The iV-copy ver- 
sion of the network model, where one expects non-trivial exponents, will be studied in a 
forthcoming publication. The scheme described in this paper can lead to a classification 
of disordered critical points that parallels the classification of sub-current-algebras Tik of 
osp(2N\2N)i. For this, the dictionary in |5(J is useful. It would be interesting to compare 



this classification with the classification based on sigma models |51|. The latter classifi- 
cation is based on discrete symmetries such as time-reversal, so it is not as strong as our 
classification of the actual critical points. 

It would be very interesting to perform more extensive numerical simulations of the 
network model that vary the relative strengths of the types of disorder and thereby see 
the phases predicted in this paper. For instance the gX + phase is characterized by the 
randomness in the flux per plaquette and the tunneling dominating over the randomness in 
the individual link phases. 

Though the models we discussed are not integrable for general anisotropic couplings, 
under the RG they can flow to isotropic current-current interactions and these are generally 
thought to be integral 



2 Note added in proof: The g — > 1/g duality of the su(2) model can be extended to the network 
model and this leads to a resolution of the flows toward the poles. This will be reported on in [33]. 



22 



VI. ACKNOWLEDGMENTS 



I especially wish to thank Denis Bernard for our collaboration on the spin Quantum Hall 
effect JHJ , and Sathya Guruswamy and Andreas Ludwig for collaboration on the disordered 
XY model |24j], both of which strongly influenced this work. I am also grateful to Phil 
Argyres, Zorawar Bassi, John Cardy, John Chalker, Vladimir Dotsenko, Bogomil Gerganov, 
Jane Kondev, Marco Moriconi, and Tony Zee for discussions. I also want to thank the 
Centre de Recherche Mathematiques at the Universite de Montreal for holding a stimulating 
workshop on related topics. This work is in part supported by the NSF. 



23 



REFERENCES 



P. W. Anderson, Phys. Rev. 112 (1958) 1900. 

E. Abrahams, P. W. Anderson, D. C. Licciardello and T. V. Ramakrishnan, Phys. Rev. 
Lett. 42 (1979) 673. 

R. B. Laughlin, Phys. Rev. B23 (1981) 5632. 

B. I. Halperin, Phys. Rev. B25 (1982) 2185. 

F. J. Wegner, Z. Phys. B35 (1979) 208; ibid B51 (1983) 579. 

A. M. M. Pruisken, Nucl. Phys. B285[FS19] (1987) 7 19; Nucl. Phys. B290[FS20] 
(1987) 61. 

P. Fendley, |cond-mat/ 0008372- 

A. A. Belavin, A. M. Polyakov and A. B. Zamolodchikov, Nucl. Phys. B241 (1984) 333. 
D. Bernard, (Perturbed) Conformal Field Theory Applied to 2D Disordered Systems: 
An Introduction. Cargese Lectures, [hep-th/9509137] . 

C. Mudry, C. Chamon and X.-G. Wen, Nucl. Phys. B466 (1996) 383. 

A. Nersesyan, A. Tsvelik and F. Wegner, Phys. Rev. Lett. 72 (1994) 2628. 
V. Dotsenko, M. Picco, and P. Pujol, Phys.Lett. B347 (1995) 113. 

A. Ludwig, Nucl. Phys. B330 (1990) 639. 

B. Gerganov, A. LeClair and M. Moriconi, On the (3 eta function for Anisotropic Current 
Interactions in 2D , |hep-th/0011189 . 



M. B. Halpern, E. Kiritsis, N. A. Obers and K. Clubok, Phys.Rept. 265 (1996) 1-138; 
|hep-th/950lT4l| . 

J. M. Kosterlitz and D. J. Thouless, J. Phys. C6 (1973) 1181; J. M. Kosterlitz, ibid C7 
(1971) 1016. 

J. T. Chalker and P. D. Coddington, J. Phys. C21 (1988) 2665. 

C.-M. Cho and J. T. Chalker, Phys. Rev. B54 (1996) 8708, |cond-mat/9605073. 



A. W. W. Ludwig, M. P. A. Fisher, R. Shankar and G. Grinstein, Phys. Rev. B50 
(1994) 7526. 

M. R. Zirnbauer, J. Math. Phys. 38, 2007 (1997), |corTcT-mat/9701024| . 

V. Novikov, M. Shifman, A. Vainshtein and V. Zakharov, Nucl. Phys. B229 (1983) 381; 

I. I. Kogan and M. Shifman, Phys.Rev.Lett. 75 (1995) 2085. 



H. H. Lin, L. Balents and M.P.A. Fisher, |cond-mat/ 980 1285 



R. M. Konik, H. Saleur and A. W. W. Ludwig, |cond-mat/ 0009 166 



S. Guruswamy, A. LeClair and A. W. W. Ludwig, Nucl.Phys. B583 (2000) 475, |concL" 
mat/ 9909143; , 



J. L. Cardy and S. Ostlund, Phys. Rev. B25 (1982) 6899; C. Mudry and X.-G. Wen, 
Nucl. Phys. B549 (1999) 613. 

Y. Hatsugai, X.-G. Wen and M. Kohmoto, Phys. Rev. B56 (1007) 1061. 

N. Hatano and D. R. Nelson, Phys. Rev. Lett. 77 (19 96) 570; Phys. Rev. B56 (1997) 

8651. 

V. Knizhnik and A. Zamolodchikov, Nucl. Phys. B247 (1984) 83. 
E. Witten, Commun. Math. Phys. 92 (1984) 455. 
D. Bernard and A. LeClair, cond-mat / 00030751 . 



P. Goddard, A. Kent and D. Olive, Commun. Math. Phys. 103 (1986) 105. 
A. B. Zamolodchikov and V. A. Fateev, Sov. Phys. J.E.T.P. 62 (1985) 215. 
D. Bernard and A. LeClair, in preparation. 



24 



D. Bernard and A. LeClair, Phys. Lett. B247 (1990) 309. 

J. Balog and A. Hegedus, J.Phys. A33 (2000) 6543, |iep-th/00032"58| . 

D. J. Amit, Y. Y. Goldschmidt and G. Grinstein, J. Phys. A13 (1980) 585. 

K. B. Efetov, Phys. Rev. Lett. 79 (1997) 491. 

D. Friedan, E. Martinec and S. Shenker, Nucl. Phys. B271 (1986) 93. 
Z. Bassi and A. LeClair, Nucl. Phys. B578 (2000) 577. 
L. Rozansky and H. Saleur, Nucl. Phys. B376 (1992) 461. 
M. Zirnbauer, |hep-th/9905054i 

V. Gurarie, Nucl. Phys. B410 (1993) 535; |hep-th/ 9303m 



H. Saleur, Nucl. Phys. B382 (1992) 486; |hep-th/91 11007 . 
H. Kausch, Curiosities at c = —2, |hep-th/95 10149) . 

M. Bershadsky, S. Zhukov and A. Vaintrob, Nucl.Phys. B559 (1999) 205, [hep 
th/990218q . 



N. Berkovits, C. Vafa and E. Witten, JHEP 9903 (1999) 018, |hep-th/99U2DU5 . 

A. LeClair, Phys.Rev.Lett. 84 (2000) 1292; |cond-mat/9905222j . 
S. Trugman, Phys. Rev. B27 (1983) 7539. 

B. Huckestein, Rev. Mod. Phys. 67 (1995) 357. 

L. Frappat, A. Sciarrino and P. Sorba, Commun. Math. Phys. 121 (1989) 457; |hcp- 
| th/9607T6T| . 

[51] A. Altland and M. Zirnbauer, Phys. Rev. B55 (1997) 1 142; M. Zirnbauer, J. Math. 
Phys. 37 (1996) 4986. 



25 



